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Abstract: We address the question of whether integrable models allow for VT- 
symmetric deformations which preserve their intgrability. For this purpose we carry 
out the Painleve test for PT-symmetric deformations of Burgers and the Korteweg-De 
Vries equation. We find that the former equation allows for infinitely many deformations 
which pass the Painleve test. For a specific deformation we prove the convergence of 
the Painleve expansion and thus establish the Painleve property for these models, which 
are therefore thought to be integrable. The Korteweg-De Vries equation does not allow 
for deformations which pass the Painleve test in complete generality, but we are able to 
qq ■ construct a defective Painleve expansion. 
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1. Introduction 



Classical as well as quantum mechanical models, which are invariant under a simultaneous 
parity transformation P : x — > —x and time reversal T : t — > —t, can be deformed in 
a controlled manner to produce new PT-symmetric theories ||, ||, ||, ||, [7|, ||. The 

crucial feature of these models is that the "PT-symmetry can be utilized to guarantee the 
reality of the energy spectra, which is due to the fact that its operator realization is a 
specific example of an anti-linear operator Q. In contrast to standard textbook wisdom, 
this means when the systems are Hamiltonian, they are non-dissipative despite being non- 
Hermitian. An important question to answer in this context is whether it is possible 
to deform models in a symmetry preserving manner. Regarding supersymmetry, it was 
recently shown |l(J that this is indeed possible. Here we will focus on the question of 
whether this is also accomplishable with regard to the symmetry underlying integrability. 
In other words, do integrable PT-symmetric models allow for deformations which do not 
destroy the integrability? A positive answer to this question will naturally lead to new 
integrable models. For some cases partial results already exist fy], (fj, ||, p], 0, 



H [nj. Here we will focus on two prototype models of integrable systems, the Burgers 
equation and the Korteweg-de Vries (KdV) equation. We will carry out the Painleve test 
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for 'PT-symmetric deformations of these models, establish thereafter in some cases the 
Painleve property and draw conclusions about their integrability. 

As there exist various notions and definitions about integrability, the Painleve test, 
the Painleve property, etc, let us briefly indicate which ones we are going to adopt in this 
manuscript. To start with, there is clearly no doubt that integrability is an extremely 
desirable property to have in a physical system, as it usually leads to exact solvability 
rather than to mere perturbative results. In the context of 1+1 dimensional quantum field 
theories the notion of integrability is usually used synonymously to the factorization of 
the scattering matrix, where the latter can be achieved simply by making use of one non- 
trivial charge [17]. Unlike as in most scenarios when one compares quantum and classical 
theories, the latter appear to be more complicated in this particular regard. In classical 
systems the definitions of integrability are much more varied and non-uniform. A common 
notion is so-called Liouville integrability, which assumes for a system with iV degrees of 
freedom the existence of N analytic single valued global integrals of motion in involution. 
The equations of motion are then separable and exact solutions can be obtained, at least 
in principle. Focussing on differential equations, as we do in this paper, one calls them 
integrable when, given a sufficient amount of initial data, they are solvable via an associated 
linear problem. The problem with all these definitions is that one does not know a priori 
whether a system is integrable or not without having computed all integrals of motion, 
mapped the problem to a linear one or actually solved the equations of motion. A general 
method to identify integrable models before this, often very difficult, task is completed 
does not exist. The closest one may get to such a method is to check whether the system 
possesses the Painleve property. One can then assume that the Painleve property implies 
integrability in the above specified sense, albeit this connection is not rigorously proven. To 
make matters worse, there exist even definitions which include the notion of integrability 
into the definition of the Painleve property [l^] . 

The concept of the Painleve property can be traced back more than a century to the 
original investigations of Painleve et al. |l!| , who set out to construct new functions from 
the solutions of ordinary differential equations (ODE). The notion of a function implies 
immediately that the solutions one is seeking ought to be single valued, which leads to a 
natural definition: An ODE whose (general) solutions have no movable 1 critical 2 singular- 
ities is said to possess the (generalized) Painleve property [18, |20|, 21]. The classification of 
possible solutions to this problem can be organised into equivalence classes obtained from 
linear fractional (Mobius) transformations and has been completed only to some degree. It 
is proven that all linear ODE posses the Painleve property, first order algebraic nonlinear 
equations lead to Weierstrass functions and second order algebraic nonlinear equations lead 
to the famous six Painleve transcendental functions. The classification of algebraic ODEs 
with Painleve property of order greater than two is still an open problem, albeit some 



partial results exist [22, p3, 24] 



The situation is somewhat less structured for partial differential equations (PDE). 
Extrapolating the previous notions one defines: A PDE whose solutions have no mov- 



1 Movable means that the solution depends on the initial values. 
2 A critical singularity is multivalued in its neighbourhood. 
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able critical singularities near any noncharacteric 3 manifold is said to possess the Painleve 
property. In general this is difficult to establish, however, there exists a more applicable 
necessary, albeit not sufficient, condition for a PDE to possess the Painleve property, which 
was developed by Weiss, Tabor and Carneval |2f| and is usually referred to as the Painleve 
test. This method is extremely practical and can be carried out in a very systematic 
fashion. Roughly speaking the main idea is that one expands the solution for a PDE (or 
ODE) in a power series starting with some single valued leading order terms. In case the 
series can be computed and involves as many free parameters as the order of the PDE then 
it is said that the PDE passes the Painleve test. In order to extrapolate from the Painleve 
test to the Painleve property one should also establish the convergence of the series, which, 
however, has been carried out only in very rare cases. 

For our purposes the relation between the Painleve property (test) and integrability 
is the most interesting. Ablowitz, Ramani and Segur (2(| conjectured almost thirty years 
ago: Any ODE which arises as a reduction of an integrable PDE, possibly accompanied by 
a variable transformation, possesses the Painleve property. To this day this conjecture has 
not been proven rigorously, but is supported by a huge amount of evidence. On one hand 
one has verified this property for almost all known integrable PDEs |2^, 27, 28, [], 29] and 



in turn, which is more impressive, one has also used it to identify new integrable ODEs 



|3C, 31]. The latter is what we hope to achieve in this manuscript. 

In summary, we will adopt here the logic that a PDE which passes the Painleve test 
and whose Painleve expansion converges also possesses the Painleve property. We take this 
as a very good indication that the system is integrable. 

We briefly explain the deformation procedure in section 2 and carry out the analysis 
for Burgers and the KdV equation in subsection 2.1 and 2.2, respectively. We state our 
conclusions in section 3. 

2. PT-symmetrically deformed integrable models 

Given a PT-symmetric PDE as a starting point, we adopt the deformation principle 
of @, H H to define new PT-symmetric extensions of this model by replacing ordinary 



derivatives by their deformed counterparts 

d x f(x) -> -i(if x ) e =: f x - e with e G R. (2.1) 

Clearly the original PT-symmetry is preserved. In general the deformations will continue 
real derivatives into the complex plane, unless e = 2n — 1 with n G 7L. We do not make use 



here of the possibility to deform also the higher derivatives via the deformation (2.1), i.e. 
replacing for instance d^f(x) by f x - £ o f x . £ , but simply define them as successive action of 
ordinary derivatives on one deformation only 

d n f{x) ^ f-lQn-l^y = d n-l fx£ = . ^ ^ 

This deformation preserves the order of the PDE. We can now employ this prescription to 

introduce new PT-symmetric models. 

3 On a characteristic manifold we can not apply Cauchy's existence theorem and therefore we do not 
have a unique solution for a given initial condition. 
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2.1 Painleve test for the "PT-symmetrically deformed Burgers' equation 

Burgers' equation is extensively studied in fluid dynamics and integrable systems, as it 
constitutes the simplest PDE involving a nonlinear as well as a dispersion term 

u t + uu x = au xx . (2.3) 

Obviously equation ( |2.3| ) remains invariant under the transformation t — ► —t, x — > —x, 
u — * u and a — > —a. Taking the constant a to be purely imaginary, i.e. a G iK, this 
invariance can be interpreted as a T^T-symmetry, which was also noted recently by Yan 
1 32]. A similar complex, albeit not PT-symmetric, version of Burgers' equations plays an 



important role in the study of two-dimensional Yang-Mills theory with an SU(N) gauge 
group 34]. The models considered in 34] become PT-symmetric after a Wick 



rotation, i.e t — > it. 

Let us now consider the "PT-symmetrically deformed Burgers' equation 

Ut + uu x - £ = iKu XX \ii with k,e,/i£1, (2-4) 

where for the time being we allow two different deformation parameters e and \i. 

Our first objective is to test whether this set of equations passes the Painleve test. 
Following the method proposed in pfifl , we therefore assume that the solution of (P-4| ) 
acquires the general form of the Painleve expansion 

oo 

u{x,t) = ^A fc (x,t)0(x,t) fc+a . (2.5) 

k=0 

Here a G Z_ is the leading order singularity in the limit (p(x, t) = (<p(x, t) — (p ) — > 0, with 
<p(x, t) being an arbitrary analytic function characterizing the singular manifold, tp being 
an arbitrary complex constant which can be utilized to move the singularity mimicking 
the initial condition and the A^(x,i) are analytic functions, which have to be computed 
recursively. 

2.1.1 Leading order terms 

As a starting point we need to determine all possible values for a by substituting the first 
term of the expansion ( |2.5D , that is u(x,t) — > Xo(x, t)<p(x, t) a , into ( |2.4[ ) and reading off 
the leading orders. For the three terms in ( |2.4| ) they are Ut ~ <p a ~ l , uu x - £ ~ anc [ 
u xx -fj, ~ (j) 01 ^^ 1 . In order for a non-trivial solution to exist the last two terms have to 
match each other in powers of (j), which immediately yields a = (e— fi— l)/(e — fi+ 1) G Z_ . 
Thus a = — 1 and e = fi is the only possible solution. This means we observe from the 
very onset of the procedure that only the models in which all x-derivatives are deformed 
with the same deformation parameter have a chance to pass the Painleve test. Therefore 
we can conclude already at this stage that one of the deformations of (^) studied in p2[ , 
i.e. e = 1 and \x generic, can not pass the Painleve test. Hence they do not possess the 
Painleve property and are therefore not integrable. 
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2.1.2 Recurrence relations 

Substituting next the Painleve expansion ( |2.5| ) for u(x,t) with a = — 1 into (^4|) with 
e = \x gives rise to the recursion relations for the A^ by identifying powers in <j)(x,t). We 
find 

at order — (2e + l): Ao + i2sK(j) x = 0, 

at order — 2e: <i>t&E,i + ^i4>x ~ i^fixx = 0> (2-6) 

at order — (2e — 1): d x ((p t 5 et i + \i4> x — iK£(j> xx ) = 0, 

such that 

Ao = — il£K(f) x , Ai = (i£K(f> xx — 4> t 5 £t i)/(j) x and A2 is arbitrary. (2-7) 

This means the number of free parameters, i.e. c/? and A2, at our disposal equals the order 
of the PDE, such that ( |2.4D passes the Painleve test provided the series ( |2.5[ ) makes sense 
and we can determine all Xj with j > 2. To compute the remaining Aj we need to isolate 
them on one side of the equation and those involving A^ with k < j on the other side. We 
expect to find some recursion relations of the form 

9(j, <I>U <Ar; <Pxx> ■ ■ = /( A J-1> A J-2, • • • , A 1; A , <j> t , <j> x , <f> xx , ■ ■ ■), (2.8) 

with g and / being some functions characteristic for the system under consideration. We 
will not present here these recursion relations for generic values of e as they are rather 
cumbersome and we shall only present the first non-trivial deformation, that is the case 
e = 2. 

2.1.3 Resonances 

For some particular values of j, say j = r%, . . . , ri, we might encounter that the function g in 
Q2.8| ) vanishes. Clearly this leads to an inconsistency and a failure of the Painleve test unless 
/ also vanishes. In case this scenario occurs, it implies that the recursion relation ( |2.£| ) 
does not fix Aj and the compatibility conditions g = f = lead to £ so-called resonances 
X n for i = 1, . . . I. When t + 1 is equal to the order of the differential equation we can in 
principle produce a general solution which allows for all possible initial values. It might 
turn out that some missing free parameters are located before the start of the expansion 
( |2.5| ), i.e. at j < 0, so-called negative resonances which can be treated following arguments 
developed in [5f|. When not enough additional free parameters exist to match the order 
of the differential equation, the series is still of Painleve type and is called defective. 

It is straightforward to determine all possible resonances by following a standard ar- 



gument. The first term in the expansion (2.5) gives rise to the leading order singularity 



which needs to be cancelled by some yet unknown term in the expansion. Let us carry out 



the calculation for Burgers equation. Using the expression for Ao from (2.7) and making 
the ansatz 

u(x,t) = -2i£K^ + # r ~ 1 , (2.9) 

4> 

we can compute all possible values of r for which •& becomes a free parameter. Substituting 
u(x,t) into ( |2.4D and reading off the terms of the highest order, i.e. _2e ~ 1+r , we find the 
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necessary condition 

i2 £ " Vtf(r + l)(r - 2)K £ 4> 2 x e = 0, (2.10) 

for a resonance to exist. This yields precisely to two resonances, one at r = 2, corresponding 
to the third equation in (j2.6|), and the so-called universal resonance at r = — 1. This means 
also at higher order we can not encounter any inconsistencies or possible breakdowns of 
the Painleve test for any value of the deformation parameter e. 

2.1.4 From the Painleve test via Painleve property to integrability 

Once it is established that a PDE passes the Painleve test one needs to be cautious about 
the conclusions one can draw as it is only a necessary but not sufficient condition for the 
Painleve property. In case one can also guarantee the convergence of the series the PDE 
possess the Painleve property, which is taken as very strong evidence for the equation to 
be integrable. This step has only been carried out rigorously in very rare cases, e.g. in 



1 36, 37]. Here we establish the convergence for one particular deformation. 



2.1.5 The e = 2 deformation 

As already mentioned, the details of the recursion relation for generic values of e are rather 
lengthy and we shall therefore only present the case e = 2 explicitly. In that case the 
deformed Burgers' equation ( |2.4| ) becomes 

ut + iuul + 2ku x u xx = (2-11) 

The substitution of the Painleve expansion ( |2.5| ) into fl2.11| ) and the subsequent matching 
of equal powers in <f> then yields the recursion relation 



i\ ct>l {Xj \{2j - 3)A - 2i{{j - 5)j + 4)k0J + 2A (A„ + 2i«0 s ) S 0J } 



(2-12) 



+ {Xj- m ~n-2Xm,xX n ;x + {m — l)X m (j) x [(n — l)Xj- m - n X n (j) x + 2Aj_ m _ n _i X n]X ] } 



n,m=l 
3-1 



-n-lX n 4> x + Xj- n -2X n ;x] ~ 2Xq(/) x [(n — \)Xj- n X n (j) x + Aj_ n _iA n;3; ] 



n=l 



2iti ^Xj— n x [A n _3;a;x 

+ (n — 3) ((n — 2)A n _i^ + 2X n -2, x 4 > x + ^n-24> xx )] 
+ (j -n - l)Xj- n (j) x [X n - 2 - xx + (n - 2) ((n - l)A n ^ + 2X n -.i :X (j) x + A n _i0 xx )] }} 
+ 2X 0:X [{j - 5)j + 6] KXj-i4>l + A i-2 [2(j - 3)k<P xx + iX 0tX ] 

~ 2Xq4> x {Xj-l [(j - 2)k4> xx + iX ; X ] + K [Xj- 2] xx + 2{j - 2)tf) x Xj-X;x]} 
+ (J - ^)Xj-3<f>t + ^J-4;i + 2kA ,x [Xj- 3 - xx + 2(j - 3)(f) x Xj-2- x ] , 

which is indeed of the general form (|2.S| ) . Having brought all Aj with j > k to the left hand 
side of (gig ), we may now successively determine the Xj to any desired order. Starting 
with the lowest value j = the equation ( 2.12j ) reduces to 



KK (Ao + i^x) = 0, 



(2-13) 
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which leads to Ao = —i^n<j) x and thus simply reproduces the expression in (|2.7D for e = 2. 
For j = 1 the equation ( 2,12j ) simplifies to 



-AqAi^I = 2X (j) x [iK\ cp xx + (A + i^n4> x ) A 0; x] , (2.14) 



such that Ai = i2n(j) xx /(j) x , which coincides with ( |2.7D for e = 2. When j = 2 the equation 
acquires the form 

A o A 2 2 (A + 4cr<£j = 2(f) x Xi. x Xl - Xl. x X + 2Xi(f) x X - x - 2iK(f) xx X - x - Uk$ x X\. x 

-2iK<p x (X 0>xx - 2cf) x Xi- x ) A . (2-15) 

It is evident that the left hand side vanishes identically and upon substitution of the values 
for Aq and Ai. We can verify that this also holds for the right hand side of (2.15), thus 



leading to the first resonance at level 2 and therefore to an arbitrary parameter A 2 . One 
may now continue in this fashion to compute the expansion to any finite order, but before 
we embark on this task we make a few further simplification. 

As the singularity has to be a noncharacteristic analytic movable singularity manifold, 
we employ the implicit function theorem and make a further assumption about the specific 
form of Xk(x, t) = Afc(i) and (fi(x, t) = x — £(£), with £(t) being an arbitrary function. Then 
the equation ( |2.12| ) simplifies to a much more transparent form 



8k 2 (8K5o,j + i(j - 2)(j + l)Xj(t)) = ^2 m)(n - l)A m (t)A i _ m _ n (t)A n (t) (2.16) 

n,m=l 

3-1 

+ [ 2 <n ~ 1) (n 2 - n - j(n - 2) + 2) X^ n (t)X n (t)] + (j - 4)A,_ 3 (i)£'(t) - A$_ 4 (t). 

n=l 

Solving this equation recursively leads to the Painleve expansion 

u{x, t) = + X 29 + |V - - ^ 4 + O^). (2.17) 



Clearly we can use ( p. 16 ) to extend this expansion to any desired order. For the ordinary 



Burgers equations, i.e. e = 1, there exist a simple choice for the free parameters, which 
terminates the expansion, such that one may generate Backhand and Cole-Hopf transfor- 



mations in a very natural way. Unfortunately ( 2.17 ) does not allow an obvious choice of 



this form. Taking for instance A 2 = yields the expansion 



AiK O 2 i£V 6 23£'£> 9 C (3 U 



10 



u{x,t) = -— + _ - — — _ + 



2 3 k 7 x 2 8 k 3 5 x 2 9 k 3 35 x 2 13 k 5 385 x 2 13 k 5 135 x 2 14 k 5 

" 2 + ifUfirW( 3 A ^ 



+ 3185 x 2 18 k 7 385 x 2 19 k 7 7 7 9 6 25 x 2 20 k 7 + 1<P j " 1 j 
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Being even more specific and assuming a travelling wave solution, the general form of the 
movable singularity is £(t) = ut, which gives 



2 3 k 7 x 2 8 k 3 35 x 2 13 k 5 3185 x 2 18 k 7 3185 x 2 21 k 9 
56Hu 6 </> 17 93^V° 6 2 5011fwV 3 



21180 25 x 2 28 k 11 3328325 x 2 32 k 13 53003575625 x 2 38 k 15 
32971wV 6 1509727^ V 9 ,„ 30 



53003575625 x 2 41 k 17 11501775910625 x 2 46 k 19 
Clearly we can carry on with this procedure to any desired order. 



+ 0(^ 3U ). (2.19) 



Convergence of the Painleve expansion Having established that the deformed Burg- 
ers equations pass the Painleve test for any value of the deformation parameter e, let us 
now see whether the obtained series converges such we may conclude that these equations 
also posses the Painleve property. It suffices to demonstrate this for some specific cases. 
Taking for this purpose A2 = 0, we can express the expansion ( 2.18| ) in the general form 

A ' 00 

u(x,t) = — ^ + 0^a n ^ n (2.20) 

^ n=l 

and employ Cauchy's root test, i.e. Y^=i In converges if and only if lim n ^oo ItJ 1 ^™ < 1, 
to establish the convergence of the series. We can easily find an upper bound for the real 
and imaginary parts of a n 

i-r, I , R e P3n-f(£ )£}£>■■■) „ n 1 o / ni\ 

\Reas n - u \ < 1 - ~ — r for u = 0,1,2, (2.21) 

I 2 3n + 4 -»T(^)\ K \ 2n - 1 V ' 

where the Pn(£') • • •) are polynomials of finite order in t, that is Yln=o ^""t"' with 

£ < 00 and to G C. The same expression holds when we the replace real part by the 
imaginary part on both sides of the inequality. We should also comment that this point of 
the proof is not entirely rigorous in the strict mathematical sense as we have only verified 
the estimate ( p. 21 ) up to order thirty. Approximating now the gamma function in ( 2.21| ) 
by Stirling's formula as n — > 00 

r (2) ~ (2.22) 

we obtain 



l/n 

lim |Rea 3 n-^ ~ — ^ — ,_ ! = 0- (2-23) 



\Rep3r, 



n ^°° ' 2 3+ ^(2vr)^e-2( 3 ^)2-^ \k\ 2 ■ 

The same argument holds for the imaginary part, such that the series ( |2.20| ) converges for 
any value of k and choices for £(t) leading to finite polynomials p n (£'> £"j £'"> • • •)■ ^ ^ s 
straightforward to repeat the same argument for A2 7^ 0. 

Alternatively we can identify the leading order term in ( 2.1l| ) and integrate the de- 
formed Burgers equation twice. In this way we change the ODE into an integral equation 

.<„«)=*, {m + f « [ i + ^ (/ W + f f • < 2 - 24 ' 
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where g(t),f(t) are some functions of integration. When discretising this equation, i.e. 
taking the left hand side to be u n +i(x,t) and replacing all the u(x,t) on the right hand 
side of this equation by u n (x,t), we may iterate ( 2.24j ) with uo(x,i) = —Ain/[x — as 



initial condition and recover precisely the expansion Q2.17). Exploiting the Banach fixed 



point theorem one may also use (2.24) as a starting point to establish the convergence of 



the iterative procedure and therefore the Painleve expansion, similarly as was carried out 
for instance in |36|, 37 1. 



Reduction from PDE to ODE Making further assumptions on the dependence of 
u(x, t) on x and t we can reduce the PDE to an ODE, and attempt to solve the resulting 
equation by integration. A common assumption is to require the solution to be of the form 
of a travelling wave u(x,t) = Q{z) = £(x — vt) with v being constant. When v is taken to 
be real, even solutions will be invariant under the original PT-symmetry. With this ansatz 
the deformed Burgers' equation for e = 2 ( p. 11 ) acquires the form 



~<z + <£ + 2<z( zz = 0. (2.25) 
When £ 2 ^ we can re-write this equation as 



dz 

which can be integrated to 



JL(c-vz+~C* + 2kC z ) = 0, (2.26) 



C(z) = e^/3(2^)V3 ^(x) + ^(x) (2 _ 2?) 
cAi{ X ) + Bi{ X ) 

with c, c being constants, x = e l7T ^(vz — c)(2vn)~ 2 ^ and Ai(x)j Bi{x) denoting Airy 
functions. 

2.2 Painleve test for the PT-symmetrically deformed KdV-equation 

The KdV-equation was found to be "PT-symmetric and was the first equation for which 
deformations have been studied 0, ^] . Next we investigate the "PT-symmetrically deformed 
version of the KdV-equation with two different deformation parameters e and fi 

ut — Quu x - £ + u XXX]ll = with e, /i G R. (2.28) 

The case fi = 1 and e generic was considered in Q and the case e = 1 and /i generic was 
studied in g. 

2.2.1 Leading order terms 

As in the previous section we substitute u(x,t) — ► Ao(x, t)cp(x, t) a into fl2.28|) in order to 
determine the leading order term. From u t ~ </> a , uu x - E ~ (f) a+a£ ~ £ an d u xxx:fl ~ Q ^~^~ 2 
we deduce a = (s — fj, — 2)/(s — fi + 1) G Z_, such that the only solution is a = — 2 with 
e = fj,. This means neither the case fi = 1 and e generic nor the case e = 1 and fi generic 
can pass the Painleve test, but the hitherto uninvestigated deformation with e = /i has at 
this point still a chance to pass it. 



-9- 



Integrable models from, VT -symmetric deformations 



2.2.2 Recurrence relations 

Substituting the Painleve expansion ( |2.5[ ) for u(x,t) with a = — 2 into ( |2.2§| ) with e = /i 
gives rise to the recursion relations for the A& by identifying powers in (f)(x, t). We compute 

order - (3e + 2): A = |e(3e + 1)$., 

order — (3e + l): Ai = — \e(2,e + 1)4> XX , 

order - 3e: A 2 = ( 4^ -3^ + 5 



2 ~ \ 0l ; + e ' 1 6^' (2.29) 



ord er - (3e - 1): A 3 = ^f^^^Lz^^j + ^ 



order - (3e - 2): A 4 = £ ^ 3e+1) ( 6 ^ ^s^eee Z x ^ ~ I ^ ^ 4>fa + fefeaZ^La 

We find that the relation at order — (3e— 2) becomes an identity only for e = 1, which makes 
us suspect that also at higher order we will not encounter compatibility conditions and 
therefore will not have enough parameters equaling the order of the differential equation. 
To test whether new compatibility conditions arise at higher levels we can use the same 
general argument as in subsection p.1.3 . 



2.2.3 Resonances 



We try once again to match the first term in the expansion (^lj) with some term of unknown 
power. Using the expression for Aq in ( |2.29| ) and making the ansatz 



u(x,t) = ie(3e + l)4 + ^ r " 2 , (2-30) 
2 <p 

we compute all possible values of r for which •& becomes a free parameter. Substituting 
u(x, t) into (2.28) and reading off the terms of the highest order, i.e. ^ _3e_2+r , we find the 



necessary condition 

e £ H) e-1 (3e + l) e ~V + 1) [6(1 + 3e) - 2(2 + 3e)r + r 2 ] ^ = 0, (2.31) 

for a resonance to exist. We observe the presence of the universal resonance at r = — 1. 
The bracket containing the quadratic term in r can be factorized as (r — r_)(r — r+) with 
r± = -(2 + 3e) ± V9e 2 - 6e - 2, such that r± G Z for 9e 2 - 6e - 2 = n 2 with n G N. For 
the solution of this equation e± = (1 ± \Jn 2 + 3)/3 to be an integer we need to solve a 
Diophantine equation 3 + n 2 = m 2 with n, m G N, which only admits n = 1 and m = 2 as 
solution. Thus the bracket only factorises in the case e = 1 into (r — 6)(r — 4). Hence, only 
in that case the system can fully pass the Painleve test. Nonetheless, we may still be able to 
obtain a defective series if all remaining coefficients Xj may be computed recursively. This 
is indeed the case as we demonstrate in detail for one particular choice of the deformation 
parameter. 

2.2.4 e = 2 deformation 



For e = fi = 2 the deformed KdV equation (2.2q) acquires the form 



ut — 6iuu x + 2iu xx + 2iu x u xxx = (2.32) 
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Since the expression become rather lengthy for generic values in the expansion we will 
present here only the case Xk(x,t) = \k(t) and 4>(x,t) = x — £(t), with being an 
arbitrary function. We find a recursion relation of the form ( [2.8] ) 



3 j-n-l 

-28i(l+j)(j 2 - 16j +42)Aj(t) = -6i£) Yl {( m ~ 2 ^ n ~ 2)A m (t)A ri (t)A i _ m _ n (t)} 

n=l m=l 

i-i 

+2i {[(7 - k)n 3 + (k - 4)kn 2 + (18 - 5k)kn + 6fc(5 + k) - 28(6 + n)]A,_ n (t)A n (i)} 

n=l 

+A;_ 6 (t) + (j-7)A;_ 5 (t). (2.33) 
The recursive solution of this equation leads to the expansion 



U(X, t) = -2 + — 37" + » - go1 , no + ' omno - oiN^ cm + 



2 156 192192 681408 73081008 216449705472 20262348288 
340915(5') V 8 + 1867RQW +0( ^ (2 M) 



23989859332927488 758331543121152 

Thus we have obtained a solution of Painleve type for the deformed KdV equation, albeit 
without enough free parameters, i.e. without the possibility to accommodate all possible 
initial values. This means we have a so-called defective series. As in the case of the deformed 
Burgers equation it is instructive to consider the series for travelling wave solutions, i.e. 
taking = ut, which yields 

. . 7 iujc^ uj 2 (f> 8 ioj 3 (f> 13 340915wV 8 391907iw 5 (/) 23 

U(X, t) = —X + — — + 1nr|1nr| + - Qno1nno - nonono r- nooonnf7/loo + 



2 156 192192 73081008 23989859332927488 56760007181706436608 



38892808841lj 



6^28 



+ 0(<p 33 ). (2.35) 



507260097462393341102260224 



Clearly we can carry on with this analysis to any desired order. The convergence of the 
expansion can be established in a similar fashion as we demonstrated for Burgers equation 
in the previous subsection or by making use of an integral equation of the type ( 2.24 ). We 
find a similar behaviour for other values of e. 



3. Conclusion 

We have carried out the Painleve test for "PT-symmetric deformations of the Burgers 
equation and the KdV equation. When deforming both terms involving space derivatives, 
we found that the deformations of the Burgers equation pass the test. In specific cases 
we have also established the convergence of the series, such that these equations have 
in addition the Painleve property. Based on the conjecture by Ablowitz, Ramani and 
Segur we take this as very strong evidence that these equations are integrable. Regarding 
these models as new integrable systems leads immediately to a sequence of interesting 
new problems related to features of integrability, which we intend to address in a future 
publication It is very likely that these systems admit soliton solutions and it should 
be possible to compute the higher charges by means of Lax pairs, Dunkl operators or 
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other methods. We should point out that most of our arguments will still hold when we 
start in ([2.4J ) with the usual Burgers equation, which has broken PT-symmetry, i.e. with 
a = in £ M. However, when embarking on the computation of charges and in particular 
energies we expect to find a severe difference as then the PT-symmetry has a bearing on 
the reality of the eigenvalues of the charges. 

For the KdV equation our findings suggest that their PT-symmetric deformations are 
not integrable, albeit they allow for the construction of a defective series. 

In future work one could also include deformations of the term involving the time 
derivative and it would clearly be very interesting to investigate other PT-symmetrically 
integrable systems in the manner in order to establish their integrability. 

Acknowledgments: P.E.G.A. is supported by a City University London research stu- 
dentship. 
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